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Abstract 

Background: The exact aetiology of canine sino-nasal aspergillosis (SNA) is unknown. In man, dysfunction in innate 
immunity, particularly in the function of pattern recognition receptors, is implicated in the pathogenesis of 
inflammatory sino-nasal disease and in fungal diseases. Associations between single nucleotide polymorphisms 
(SNPs) in Toll-like receptors (TLRs) and these diseases have been identified. Similarly, in dogs SNPs in genes 
encoding TLRs may be important in the pathogenesis of SNA. The aims of the present study were (1 ) to identify the 
presence of non-synonymous SNPs in the coding regions of the TLR2, 4 and 9 genes in dogs suffering from SNA, and 
(2) to investigate the SNP genotypes in dogs with SNA compared with a control population. 

Results: Direct sequencing of nine dogs of various breeds with SNA revealed two non-synonymous SNPs in the coding 
region of TLR2, eight in TLR4 and four in TLR9. These non-synonymous SNPs were further evaluated in a case-control 
study of affected Golden Retrievers, Labrador Retrievers, Rottweilers and Beaucerons. Genotyping was performed using a 
combination of allele-specific primers and hydrolysis probe assays in 31 dogs with SNA and 31 controls. No significant 
difference in minor allele frequency was identified between these groups, for all studied SNPs, in any of the four breeds. 

Conclusions: These findings do not support a role for non-synonymous SNPs in the TLR 2, 4 and 9 coding regions in the 
pathogenesis of canine SNA, but do not exclude a role for innate immunity in the pathogenesis of the disease. 
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Background 

The sino-nasal mucosa is an important interface with 
the environment and represents a dynamic system for 
innate host defence [1]. In addition to providing a bar- 
rier to the entry of pathogens, epithelial cells lining the 
airway are able to recognise pathogen-associated mo- 
lecular patterns (PAMPs) [2]. Toll-like receptors (TLRs) 
are situated on these cells and on dendritic cell popula- 
tions in this location and appear to play crucial roles in 
the distinction between commensal flora and PAMPs. 
Single nucleotide polymorphisms (SNPs) within the TLR 
genes have been reported to impact the immune reac- 
tion to microbial antigens, and to be associated with the 
incidence and outcome of infections [3]. Collectively, 
the TLRs cover the recognition of a wide variety of 
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pathogens (i.e. viruses, bacteria, parasites and fungi). 
TLR2, TLR4 and TLR9 are particularly associated with rec- 
ognition of fungi (e.g. Candida albicans and Aspergillus 
fumigatus) [4-9], and SNPs in these TLRs have been asso- 
ciated with higher risk of fungal disease in man [10]. 

In dogs, sino-nasal aspergillosis (SNA) is a cause of 
chronic rhinitis, and the disease is most often caused by 
A. fumigatus [11]. SNA is characterized by the formation 
of superficial mucosal fungal plaques within the nasal 
cavity and/or frontal sinus, associated with a severe mu- 
cosal inflammatory response and a local destruction of 
nasal bone [12]. As the disease usually affects systemic- 
ally healthy dogs, and A. fumigatus is a ubiquitous fun- 
gus, a local immune deficiency in affected dogs has been 
suggested [13]. 

The clinicopathological features of SNA are similar to 
those of human chronic erosive non-invasive fungal si- 
nusitis [14]. Although TLRs have not been studied in 
this type of sinusitis in people, they have been better 
characterized in chronic rhinosinusitis (CRS) including 
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fungal CRS, and some forms of aspergillosis. Pathogen- 
esis of the CRS complex remains unclear, but association 
of this disease with TLR polymorphisms, particularly 
TLR2 SNPs, has been described [15]. 

Association between TLR SNPs and incidence of vari- 
ous forms of human Aspergillus infection have been re- 
ported. A significant association between the incidence 
of invasive aspergillosis (IA) in recipients of transplanted 
allogeneic stem cells and the presence of some TLR4 
SNPs in donors has been reported [16,17]. Hence, 
Carvalho et al. [10] have described an increased fre- 
quency of chronic pulmonary aspergillosis among TLR4 
SNP carriers. In the same study, a SNP in TLR9 was 
shown to predispose to the allergic form of pulmonary 
aspergillosis. 

Polymorphisms in TLR genes have not yet been de- 
scribed in dogs with SNA, but alteration in these recep- 
tors could hamper their function, therefore disrupting 
the recognition of fungi, eventually leading to fungal es- 
cape from immune surveillance. Recently, up-regulation 
of TLR2, 4 and 9 mRNA expression has been reported 
in the nasal mucosa of affected dogs [18]. Comparison 
of the canine and human TLR genetic sequences show a 
high degree of homology [19]; suggesting the polymor- 
phisms in TLR genes found in people with aspergillosis 
and/or fungal CRS might also be observed in dogs. 

A retrospective analysis of the medical records at 
the Veterinary Clinic of the University (CVU) of Liege 
(Belgium) over the last 5 years revealed four predominant 
canine breeds diagnosed with SNA: Golden Retriever, 
Labrador Retriever, Rottweiler and Beauceron Sheepdog. 
A candidate gene approach was used to test the hypoth- 
esis that SNPs in TLR2, 4 and 9 genes were associated 
with canine SNA. The presence of non-synonymous 
SNPs in the coding regions of the TLR2, 4 and 9 genes 
was investigated by direct sequencing of PCR-amplified 
products from these regions, using DNA from nine dogs 
with SNA belonging to various breeds (including the 
four aforementioned breeds). Secondly, a combination of 
allele-specific primers and Taq-man probe-based PCR 
genotyping was applied to samples from affected and 
control dogs derived from these four breeds in order to 
assess the significance of these SNPs in canine SNA. 

Results 

Selection of single nucleotide polymorphisms 

The initial part of the study focused on a mutational 
analysis of TLR2, 4 and 9 genes in nine dogs with SNA 
belonging to various breeds. This analysis revealed two 
non-synonymous SNPs in the coding region of TLR2: 
C137A (rs22410121) and C1547T (not reported in the 
canine genome database). Eight non-synonymous SNPs 
were identified in the TLR4 exons: T23C (not reported in 
the canine genome database but described by Kathrani 



et al. [20]), C107T (not reported previously), TS00C 
(rs22145736), A600C {rs22189454), A688G (rs221894S6), 
G 1039 A (not reported in the canine genome database but 
described by Kathrani et al. [20]), A1S71T (rs22124023) 
and G1807A (rs22123995). Four non-synonymous SNPs 
were identified in TLR9 exon: G1138A (not reported 
previously), A1372C (not reported previously), A2158G 
\rs22882109) and G2927A (not reported previously). 

All non-synonymous SNPs identified resulted in a 
change in the class of amino acid coded (Table 1). 

SNP genotyping and association study 

To assess the significance of the identified SNPs, the allelic 
frequencies in a group of dogs with SNA belonging to four 
breeds commonly affected by SNA (Golden Retriever, 
Labrador Retriever, Rottweiler and Beauceron Sheepdog) 
were compared with those found in control dogs of the 
same breeds (Table 2), and no significant difference was 
found. Two SNPs (TLR2 C1547T and TLR4 TS00Q were 
found monomorphic between and within samples from 
the four selected breeds. 

The same analysis was conducted within the four 
breed groups separately. Again, the allelic differences, 
although variable across breeds, were not related to 
the disease phenotype (Tables 3, 4, 5 and 6). In Golden 
Retrievers, the TLR4 C107T SNP was not found, and the 
least common alleles for all breeds taken together for 
TLR9 A1372C and TLR9 A2158G were the most common 
in this breed. In Golden Retrievers, the association of the 
C137A SNP with the disease almost reached significance. 
In the Labrador Retriever dogs, the least common alleles 

Table 1 Amino acid change coded by non-synonymous 
single nucleotide polymorphisms in the canine TLR2, 4 
and 9 genes 



Gene SNP (accession number) Amino acid Amino acid change 







wild type 


associated with SNP 


TLR2 


C137A (rs22410121) 


Serine 


Tyrosine 




C1547T 


Serine 


Leucine 


TLR4 


T23C 


Valine 


Alanine 




C107T 


Threonine 


Isoleucine 




T500C (rs22 145736) 


Leucine 


Proline 




A600C (rs22 189454) 


Histidine 


Glutamine 




A688G (rs22 189456) 


Lysine 


Glutamic acid 




G1039A 


Alanine 


Threonine 




A157lT(rs22124023) 


Glutamic acid 


Valine 




G1807A [rs221 23995) 


Glutamic acid 


Lysine 


TLR9 


G1138A 


Glutamic acid 


Lysine 




A1372C 


Threonine 


Proline 




A2158G {rs22882109) 


Serine 


Glycine 




G2927A 


Arginine 


Histidine 
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Table 2 Single nucleotide polymorphism allele association 
with sino-nasal aspergillosis dogs 



SNP Minor Minor allele Minor allele p value 

allele frequency in the frequency in the 
31 SNA dogs 31 CTL dogs 



TLR2 C137A 


C 


0.42 


0.48 


0.47 


TLR2 C1547T 


T 


0 


0 


NA 


TLR4 T23C 


T 


0.40 


0.40 


1 


TLR4 C107T 


T 


0.03 


0 


0.15 


TLR4 T500C 


C 


0 


0 


NA 


TLR4 A600C 


C 


0.48 


0.43 


0.59 


TLR4 A688G 


A 


0.47 


0.43 


0.71 


TLR4 G1039A 


A 


0.34 


0.34 


I 


TLR4A1571T 


A 


0.13 


0.10 


0.5/ 


TLR4 G1807A 


G 


0.13 


0.10 


0.57 


TLR9G1138A 


A 


0.06 


0.05 


0.70 


TLR9A1372C 


A 


0.50 


0.37 


0.15 


TLR9 A2158G 


G 


0.55 


0.43 


0.21 


TLR9 G2927A 


A 


0.03 


0 


0.15 



Statistical significance was set at P < 0.05. 



for all breeds taken together for TLR2 C137A, TLR4 
A600C, TLR4 A688G, TLR4 and G1039A were most com- 
mon. In Rottweiler dogs, TLR4 C107T, TLR4 A1571T, 
TLR4 G1807A, TLR9 G1138A and TLR9 G2927A SNPs 
were absent, and the least common allele for all breeds 
taken together for TLR4 A600C was the most common 



Table 4 Single nucleotide polymorphism allele association 
with sino-nasal aspergillosis Labrador Retriever dogs 



SNP 


Minor 
allele 


Minor allele 
frequency in the 
SNA Labrador 
Retriever dogs 


Minor allele 
frequency in the 
CTL Labrador 
Retriever dogs 


p value 


TLR2 Cm A 


A 


0.50 


0.37 


0.72 


TLR2 C1547T 


T 


0 


0 


1 


TLR4 T23C 


T 


0.44 


0.25 


0.46 


TLR4 C107T 


T 


0.06 


0 


1 


TLR4 T500C 


C 


0 


0 


1 


TLR4 A600C 


A 


0.37 


0.25 


0.70 


TLR4 A688G 


G 


0.37 


0.25 


0.70 


TLR4 G1039A 


G 


0.50 


0.31 


0.47 


TLR4A1571T 


A 


0.12 


0.06 


1 


TLR4 G1807A 


G 


0.12 


0.06 


1 


TLR9G1138A 


A 


0.18 


0.06 


0.60 


TLR9A1372C 


A 


0.37 


0.19 


0.43 


TLR9 A2158G 


G 


0.56 


0.37 


0.48 


TLR9 G2927A 


A 


0.12 


0 


0.48 



Statistical significance was set at P < 0.05. 



allele. In the Beauceron Sheepdogs, TLR2 C137A, TLR9 
G1138A, TLR9 A2158G and TLR9 G2927A SNPs were 
absent, and the least common allele in all breeds taken 
together for TLR9 A1372C was the most common allele in 
this breed. 



Table 3 Single nucleotide polymorphism allele association Table 5 Single nucleotide polymorphism allele association 
with sino-nasal aspergillosis Golden Retriever dogs with sino-nasal aspergillosis Rottweiler dogs 



SNP 


Minor 
allele 


Minor allele 
frequency in 
the SNA Golden 
Retriever dogs 


Minor allele 
frequency in 
the CTL Golden 
Retriever dogs 


p value 


SNP 


Minor 
allele 


Minor allele 
frequency 
in the SNA 
Rottweiler dogs 


Minor allele 
frequency 
in the CTL 
Rottweiler dogs 


p value 


TLR2 C137A 


C 


0.15 


0.45 


0.08 


TLR2 C137A 


C 


0.39 


0.17 


0.26 


TLR2 C1547T 


T 


0 


0 




TLR2 C1547T 


T 


0 


0 


1 


TLR4 T23C 


T 


0.35 


0.40 




TLR4 T23C 


C 


0.44 


0.44 


1 


TLR4 C107T 


T 


0 


0 




TLR4 C107T 


T 


0 


0 


1 


TLR4 T500C 


C 


0 


0 




TLR4 T500C 


C 


0 


0 


1 


TLR4 A600C 


C 


0.30 


0.25 




TLR4 A600C 


A 


0.39 


0.56 


0.51 


TLR4 A688G 


A 


0.30 


0.25 


0.10 


TLR4 A688G 


A 


0.56 


0.45 


0.74 


TLR4 G1039A 


A 


0.05 


0.05 




TLR4 G1039A 


A 


0.56 


0.45 


0.74 


TLR4A1571T 


A 


0.25 


0.20 




TLR4A1571T 


A 


0 


0 




TLR4 G1807A 


G 


0.25 


0.20 




TLR4 G1807A 


G 


0 


0 




TLR9G1138A 


A 


0.05 


0.10 




TLR9G1138A 


A 


0 


0 




TLR9A1372C 


C 


0.10 


0.25 


0.41 


TLR9A1372C 


A 


0.06 


0.06 




TLR9 A2158G 


A 


0.20 


0.25 




TLR9A2158G 


G 


0.06 


0.06 




TLR9 G2927A 


A 


0 


0 




TLR9 G2927A 


A 


0 


0 





Statistical significance was set at P < 0.05. 



Statistical significance was set at P < 0.05. 
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Table 6 Single nucleotide polymorphism allele association 
with sino-nasal aspergillosis Beauceron sheepdogs 

SNP Minor Minor allele Minor allele p value 

allele frequency in the frequency in the 



SNA Beauceron CTL Beauceron 
Sheepdogs Sheepdogs 


TLR2 C137A 


A 


0 


0 


1 


TLR2 C1547T 


T 


0 


0 


1 


TLR4 T23C 


T 


0.12 


0.37 


0.57 


TLR4 C107T 


T 


0.12 


0 




TLR4 T500C 


C 


0 


0 




TLR4 A600C 


C 


0.37 


0.25 




TLR4 A688G 


A 


0.37 


0.25 




TLR4 G1039A 


A 


0.25 


0.12 




TLR4A1571T 


A 


0.12 


0.12 




TLR4 G1807A 


G 


0.12 


0.12 




TLR9G1138A 


A 


0 


0 




TLR9A1372C 


C 


0.25 


0.50 


0.61 


TLR9 A2158G 


A 


0 


0.37 


0.20 


TLR9 G2927A 


A 


0 


0 


1 


Statistical significance was set at P < 0.05. 



Discussion 

The present study investigated the association between 
SNPs in TLR2, 4 and 9 genes, and the risk of developing 
SNA in dogs. In the first part of the study, sequencing of 
the coding regions of these genes in nine dogs with SNA 
(from various breeds) revealed two, eight and four non- 
synonymous SNPs in the coding region of TLR2, TLR4, 
and TLR9 genes, respectively. In the second part of the 
study, no association was found between these SNPs and 
the presence of SNA. 

Some of the SNPs identified in the first part of the 
study were not found in samples from the four breeds of 
dog genotyped in the second part, and minor allele fre- 
quencies varied extensively between breeds, with some 
uncommon alleles in one breed becoming the most com- 
mon alleles in the other, and vice versa. This supports a 
wide genetic variations between canine breeds [21]. 

The only result close to significance was found for 
TLR2 C137A SNP in the Golden Retrievers. Due to the 
unpredictable direction of allelic frequency shift between 
dogs with SNA and controls, PLINK software was used 
to compute two-tailed Fisher's exact test p-vahies. Using 
a one tailed test would have led to half the actual jj-value 
(0.04; data not shown), thereby reaching significance. Al- 
though this might be interpreted as a significant result, 
the rationale for choosing a one-tailed test might be 
debated. Furthermore, multiple testing issues (i.e. per- 
forming several tests) could be considered, leading to a 
lowered significance threshold that would make this re- 
sult not significant. On the other hand, the power of the 



test might be low due to the relatively small size of the 
sample. Therefore, this SNP could merit further investi- 
gation with a larger number of dogs. 

Although some researchers have studied the associ- 
ation between TLR SNPs in human CRS and some fungal 
diseases, the results are heterogeneous and sometimes 
contradictory. Several researchers have shown a significant 
association between TLR4 and/or TLR9 SNPs and the inci- 
dence of several forms of aspergillosis [10,16,17]; however, 
another study failed to associate these polymorphisms with 
that disease. Recently, Carvalho et al. noted a positive asso- 
ciation between TLR4 SNPs and fungal colonization, but 
not susceptibility to fungal infection [22]. Kesh et al. did 
not report any positive association of recipient or donor 
TLR4 SNPs with the incidence of IA in patients undergo- 
ing allogeneic stem cell transplantation [23]. In human 
CRS, there are also discrepancies between studies regard- 
ing the role of SNPs in TLR genes. Park et al. reported an 
association between TLR2 SNPs in a Korean population 
[15], but Tewfik et al. could not confirm these results in 
Quebec [24]. 

The effect of the environment on susceptibility to cer- 
tain diseases is another important aspect that should be 
considered together with genetic background. Vercelli 
[25] reported that genes and environment are inter- 
twined in complex nonlinear relationships. Therefore, 
the same genetic background may result in the expres- 
sion of different phenotypes when exposed to different 
environments. Consequently, in a disease with a com- 
plex pathogenesis, the genotype of multiple contributing 
genes, together with environmental factors, may modify 
the likelihood of an individual developing the disease. In 
such instances, individual gene mutations are insufficient 
to cause the disease on their own [26]. 

There are several limitations to the present study. 
Firstly, although we cannot confirm an association be- 
tween polymorphism in the TLR2, 4 and 9 genes with 
SNA, it is possible that because of the small size of the 
study population, a slight/ weak but real contribution of 
TLR2, 4 or 9 SNPs has been missed. Another potential 
limitation is the selection of control dogs. Although ef- 
forts were made to exclude dogs with respiratory, in- 
flammatory, immune-mediated, neoplastic or infectious 
diseases, the possibility cannot be excluded that these dogs 
might develop such disease later in life. Nevertheless, to try 
to overcome this problem, we recruited only dogs older 
than 6 years. Another limitation was the use of a candidate 
gene approach, which permits comparison of allele fre- 
quencies of SNPs in genes suspected to be involved in 
SNA among dogs with the disease and others without 
(controls), but this method precludes finding other genes 
of interest. Genome-wide approaches would overcome this 
limitation, but the need for a large number of affected dogs 
and the high associated expense renders this difficult. 
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Conclusions 

The findings of the present study suggest that TLR2, 4 
and 9 are not implicated in the pathogenesis of SNA in 
dogs. However, these results are limited by the size of 
the study population. This does not rule out the possibil- 
ity of innate immune defects contributing to the devel- 
opment of the disease. Other TLR genes and/or other 
genes involved in mucosal immunity may be important 
and require further investigation. 

Methods 

Selection of single nucleotide polymorphisms 

A mutational analysis of the TLR2, 4 and 9 exons was 
carried out in nine dogs diagnosed with SNA in order to 
determine the presence of non-synonymous SNPs in the 
coding regions of the genes and to investigate these fur- 
ther in a case-control association study to determine 
their significance in SNA. 

The nine selected dogs (aged 10 months to 14 years; 
median 6.7 years; 3 females and 6 males) were from vari- 
ous breeds: Golden Retriever (n = 2), Labrador Retriever 
(n = 2), Rottweiler (n = 2), and one each of Beauceron, 
German Shepherd Dog and crossbreed. Diagnosis of 
SNA was made on the basis of physical examination, 
observation of intranasal fungal plaques associated with 
severe turbinate destruction on rhinoscopy, presence of 
serum antibody specific for A. fumigatus and culture of 
A. fumigatus from suspected lesions. A computed tom- 
ography examination of the head was performed in some 
dogs. 

Blood stored in ethylenediaminetetraacetic acid (EDTA) 
or biopsy tissue from the nasal mucosa were available for 
each dog and were used for extraction of genomic DNA 
using the Macherey-Nagel NucleoSpin blood (or tissue) 
isolation kit (ABgene, Epsom, UK). All biological material 
used to perform the present study was collected for diag- 
nostic purpose in all dogs. Nevertheless, the right to use it 
was obtained by owner consent. 

Canine TLR2, 4 and 9 forward and reverse primers 
were designed via Primer 3 (http://bioinfo.ut.ee/primer3- 
0.4.0/primer3/), using the canine specific GenBank se- 
quences, as described previously [27] (Table 7). Only 
primers for exon 3 were design for TLR9, because only two 
bases of exon 2 are coding. Primers were synthesized by 



Eurogentec SA. (Seraing, Belgium) and reconstituted in 
TE buffer (10 mM Tris, 1 mM EDTA, pH 8) before use. 

Conventional polymerase chain reaction (PCR) was per- 
formed using Taq Platinum 8 PCR master mix (Invitrogen, 
California, USA), with reactions heated to 94°C for 2 min, 
followed by 35 cycles of 94°C for 25 sec, 60°C for 25 sec 
and 68°C for 4 min then a further 72°C for 7 min. The 
PCR products obtained were separated by 2% w/v agarose 
gel electrophoresis and appropriately sized products 
were excised with a clean scalpel blade and purified 
using the Nucleospin 8 Extract II (Macherey-Nagel, Hoerdt, 
France) as per the manufacturer's instructions. Gel-purified 
PCR products were sent to GIGA research centre in Liege, 
Belgium for sequencing using sequencing primers specific- 
ally designed to cover the entire sequence of the exons 
(Table 8). 

The sequence data were compared with the canine- 
specific GenBank sequences using GeneStudio contig 
assembly editor (http://www.genestudio.com/), and non- 
synonymous SNPs were identified. 

SNP genotyping 

The non-synonymous SNPs identified in the first part of 
the study were further analysed in case and control pop- 
ulations by genotyping. 

The case group consisted of 31 dogs diagnosed with 
SNA (aged 10 months to 12 years; median 6.5 years; 
7 females and 24 males) from four different breeds: 
Golden Retriever (n = 10; minimum age 10 months; 
maximum age 12 years; median 7 years; 1 female and 
9 males), Labrador Retriever (n = 8; minimum age 1.3 years; 
maximum age 12 years; median 7.1 years; 2 females and 
6 males), Rottweiler (n = 9; minimum age 3 years; max- 
imum age 9 years; median 5 years; 2 females and 7 males) 
and Beauceron (n = 4; minimum age 10 months; maximum 
age 10.5 years; median 2.5 years; 2 females and 2 males). 
Diagnosis of SNA was based on the criteria given above. 
These dogs were presented at the Veterinary Clinic of the 
University (CVU) of Liege (Belgium) between July 2007 
and June 2013. During that period, a total of 72 dogs from 
various breeds were diagnosed with SNA at the same insti- 
tution. The breeds of dogs included in the present study 
were the 4 most frequently encountered breeds in these 72 
dogs. A Fisher's exact test was used to assess whether those 



Table 7 Primers used for amplifying 71/72, 4 and 9 for sequencing 



Gene Accession number Product length Forward primer (5'-3') Reverse primer (5'-3') 



TLR2 


NC_0065973 




2538 


AATAAACTGTCAAAACAATCACTCA 


AAGCAAGTCTGCAAAGGACA 


TLR4 


NC_0065933 


exon 1 


359 


ACAAAAGCCCAGAACGCTAA 


TGCACAGAGAGCAG I I I I ICA 






exon 2 


397 


GAGAGAGGGCAGTTGAGGTG 


AGATGAGGCAATGGGATCTG 






exon 3 


3111 


CTGAATCTGTGGGGCTTCrr 


GAGACATGAAAAATGAGAACTGGA 


TLR9 


NC_0066023 


exon 3 


3390 


GGAACCCTGTTGGGAGACC 


TTTGGGAAGGAAAGCCTGAC 
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Table 8 Additional internal sequencing primers used 

Gene Primer set Primer sequence {5'-3') 

TLR2 TLR2 sequencing forward IGGTrCCTTGCTCACTrTCA 

TLR2 sequencing reverse TGGCAAAATCAGGGAAAATG 

TLR4 TLR4 exon 3 sequencing forward GGTGTGCCAGGAATCATTTG 

TLR4 exon 3 sequencing reverse TAGGATCTGGAGGGAGAGGAG 

TLR9 TLR9 exon 3 sequencing forward GTTCAGCCGGAGATGTTTGT 

TLR9 exon 3 sequencing reverse CGAGCTTGTTATGGGAGAGG 



breeds were predisposed to the disease as compared with 
the general population of dogs presented to the CVU dur- 
ing the same period of time. The four canine breeds were 
proved to be predisposed to the disease (p < 0,01 in the 
Golden Retriever and the Rottweiler; p < 0,05 in the 
Labrador Retriever and the Beauceron). 

The control group consisted of 31 dogs of the same 
breeds: Golden Retriever (n = 10; minimum age 6 years; 
maximum age 14 years; median 9 years; 4 females and 
6 males), Labrador Retriever (n = 8; minimum age 9 years; 
maximum age 14 years; median 12 years; 3 females and 
5 males), Rottweiler (n = 9; minimum age 6 years; max- 
imum age 8 years; median 7.8 years; 4 females and 5 
males) and Beauceron (n = 4; minimum age 8 years; max- 
imum age 12 years; median 11 years; 1 female and 3 
males), with no history or signs of respiratory disease. In- 
flammatory, immune-mediated, neoplastic or infectious 
diseases were also excluded. Twelve dogs were healthy; 
five were presented for laryngeal paralysis, four for de- 
generative osteoarthritis, three for a lipoma and two for 
perineal hernia. One dog each had vaginal haemorrhage, 
arrhythmia, anal sac impaction, hyperadrenocorticism or 
hypoadrenocorticism. All dogs in this group were older 
than 6 years of age (6 - 14 years; median 9 years; 15 fe- 
males and 16 males). Owner consent was again obtained 
for each dog. 

Blood stored in ethylenediaminetetraacetic acid (EDTA) 
was used for genomic DNA extraction using the 
Macherey-Nagel NucleoSpin Blood Isolation Kit (ABgene, 
Epsom, UK). 

The allele-specific primer-based PCR was used to 
genotype the SNPs [28,29]. Briefly, allele-specific PCR 
primers were designed so that the 3' terminal nucleo- 
tide of a primer corresponded to the site of the SNPs 
(Table 9). An allele-specific primer matches perfectly 
with one allele (the specific allele) or has a 3 ' mis- 
match with a non-specific allele. Because mismatched 
3' termini are extended by Taq polymerases with much 
lower efficiency than correctly matched termini, the 
allele-specific primer preferentially amplified the specific 
allele. To improve the specificity of traditional allele- 
specific primers, we followed a modified method [30,31] 
in which an artificial base pair mismatch was introduced 



between the third and the fifth base from the 3 ' terminus 
to each primer (Table 9). Two positive controls were de- 
signed for each SNP, to include allele-specific target se- 
quence (without the addition mismatch) and the common 
primer (forward or reverse as appropriate), to test the abil- 
ity of the assay to identify the SNP (Table 9). Primers and 
positive controls were synthesized by Eurogentec S. A. 
(Seraing, Belgium) and reconstituted in TE buffer (10 mM 
Tris, 1 mM EDTA, pH 8) before use. 

The allele-specific reactions were performed by quan- 
titative PCR (qPCR) using SYBR Green I with separate 
reactions for each primer set (two sets per SNP). SNP 
genotyping was performed by comparison of Ct values 
(Figure 1) with a homozygous sample having a lower Ct 
with the specific primer set and a higher Ct with the 
non-specific one, and a heterozygous sample having ap- 
proximately equal results from both assays. Each assay 
was optimized with the synthetic positive controls and 
residual DNA from the samples used in the sequencing 
study. The qPCR was performed using GoTaq Colourless 
Master Mix (Promega Corporation, Southampton, UK) 
in a Stratagene Mx3005P (Agilent Technologies, Stockport, 
Cheshire, UK) with results analysed using the MxPro 
Software (ver. 4.1; Agilent Technologies). Optimisa- 
tion involved varying the annealing temperature (55, 
60 and 65°C), with and without additional MgCl 2 (final 
concentration: 3.0 and 4.5 mM), to balance the Ct value 
difference between the target and non-target allele-specific 
PCR, and non-specific product formation with genomic 
DNA and water negative controls. The final reactions con- 
sisted of 200nM of each primer, 4.5 mM Mg and 1:50,000 
SYBR Green I (Invitrogen) with a thermocycling protocol 
of incubation at 95°C for 2 min for Taq activation and 
then 35 cycles of 95°C for 15 sec and annealing (55 or 
60°C: Table 9) for 30 sec during which the fluorescence 
data were collected. The thermocycling protocol was ex- 
tended by heating samples from 75°C to 95°C in 0.5°C 
increments for 15 sec, during which time the fluores- 
cence data were collected in order to create a melt curve. 
For each SNP, the two genotyping assays were run for a 
single DNA sample on the same plate which also in- 
cluded a number of additional controls including the 
target and non-target synthetic positive control, and AE 
buffer (negative control). 

For three SNPs (TLR9 A1372C, TLR9 C2607T, TLR9 
G2927A) it was not possible to genotype on the basis of 
Ct value alone because Ct values were too similar with 
both homozygous and heterozygous animals due to 
the presence of additional amplified products in both 
reactions (Figure 2). Therefore, capillary electrophoresis 
(QIAxcel Advanced System, QIAGEN Group, UK) using 
the QIAxcel DNA High Resolution Kit (QIAGEN), QX 
Alignment Marker 15 bp/500 bp (QIAGEN) and QX 
DNA Size Marker 25-500 bp (QIAGEN) was used to 
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Table 9 Allele-specific primers and positive control used for genotyping 



SNP 



Primer Primer sequence (5'-3') 



Positive control 



TLR2 

C137A 

(110 bp) [60°C] 

C1S47T 

(81 bp) [60°C] 

TLR4 

C107T 

(101 bp) [60°C] 
T500C 

(58 bp) [60°C] 
A600C 

(110 bp) [60°C] 
A688G 

(97 bp) [60°C] 

G1039A 

(66 bp) [60°C] 

A1571T 

(163 bp) [55°C] 

G1807A 

(53 bp) [60°C] 

TLR9 

G1138A 

(152 bp) [55°C] 

A1372C 

(95 bp) [60°C] 

A2158G 

(58 bp) [60°C] 

G2927A 

(113 bp) [55°C] 



Forward CCGCTCCAGATCTTTGAAgTC 
CCGCTCCAGATCTTTGAAgTA 
Reverse AATCCCGAAGGTCACTGTTG 
Forward GCCTCCTTCTTACCCACCTT 
Reverse GGAATCCAGTFGCTCCTaCG 
GGAATCCAGTFGCTCCTaCA 

Forward GGAAAGGAGAGAGGGCAGTT 
Reverse AGCTCCATGCATFGGTAAGaAG 

AGCTCCATGCATFGGTAAGaAA 
Forward CAATCTTATCCATTCCTTCAAcCT 

CAATCTTATCCATTCCTTCAAcCC 
Reverse CTCCAGGTTGGGCATGTTAG 
Forward CTAACATGCCCAACCTGGAG 
Reverse TAAAGAAAGGTTGAGTAGTGGQTT 

TAAAGAAAGGTTGAGTAGTGGQTG 
Forward GGTFCCTTTAAAGAAATTAAACTCgATA 

GGTFCCTTTAAAGAAATrAAACTCgATG 
Reverse AGCCAGCCAGACCTFGAATA 
Forward ATGGCAACGGTTGGAAATAG 
Reverse AGAGAGTCCAGCTCCCAaGC 

AGAGAGTGCAGCTCCCAaGT 
Forward ATGGCTGATAATTCCTTTCC 
Reverse GGCTATTGTGACTCATATTTAaCA 

GGCTATTGTGACTCATATTTtTCT 
Forward TGGGTCAAAGACCACAGACA 
Reverse GCACACACCATTTGTTCAAgTTC 

GCACACACCATTTGTTCAAgm 

Forward GCGCAGACTCAACCTGTCCTT 
Reverse GCGACTGGAGCGTGGTgTC 

GCGACTGGAGCGTGGTgTT 
Forward GCGGCTGCCACAGGGGAGGT 
Reverse GAGCCGGGGGTGCCCAGaGT 

AGCCGGGGGTGCCCAGaGG 
Forward CTGGACCTCAGCGGCAACAG 
Reverse CGGCCAGGGCAAAAAAcCT 

CGGCCAGGGCAAAAAAcCC 
Forward CTCCTGCGTGCCAGCTTCCT 
Reverse AGCCGCACATAGCGGGtGC 

AGCCGCACATAGCGGGtGT 



CCGCTCCAGATCTTTGAAGTCGATGCATTGGCAACAGTGACCTTGGGGATT 
CCGCTCCAGATCTTTGAACTACATGCATTGGCAACAGTGACCITCGGGATT 



GGAATCCAGTTGCTCCTTCGAGAAAGAATTTGTAAGGTGGGTAAGAAGGAGGC 
GGAATCCAGTTGCTCCTTCAAGAAAGAATTTGTAAGGTGGGTAAGAAGGAGGC 



AGCTCCATGCATTGGTAAGTAGGCATGCACCTCAACTGCCCTCTCTCCTTTCC 
AGCTCCATGCATTGGTAAGTAAGCATGCACCTCAACTGCCCTCTCTCCnTCC 
CAATGTATCCATrCCTTCAAGCTTGCGTATTTCTCTAACATGCCCAACCTGGAG 
CAATCTTATCCATTCCrrCAAGCCTGCGTATrTCTCTAACATGCCCAACCTGGAG 



TAAAGAAAGGTTGAGTAGTGGCATTTGATGTAGACGTTCTCCAGGTTGGGCATGTTAG 
TAAAGAAAGGTTGAGTAGTGGCATGTGATGTAGACGTrCTCCAGGTTGGGCATGTTAG 
GGTrCCnTAAAGAAATFAAACTCCATAAACTGACTlTGAGAAGCCAGCCAGACCTTGAATA 
GGTTCCTTTAAAGAAATTAAACTCCATGAACTGACTTTGAGAAGCCAGCCAGACCTTGAATA 



AGAGAGTCGAGCTCCGATGCCTAAGTTACAATTAACTATTTCCAACCGTTGCCAT 
AGAGAGTCCAGCTCCCATGTCTAAGTTACAATTAACTATTTCCAACCG7TGCCAT 

GGCTATTGTGACTCATAT1TATCAAGGGAGTTGTCCGGAAAGGAATTATCAGCCAT 
GGCTATTGTGACTCATATTTATCTAGGGAGTTGTCCGGAAAGGAA7TATCAGCCAT 

GCACACACCATTTGTTCAACTTCCACCAAGAGCTGTCTGTGGTCTTTGACCCA 
GCACACACCATTTGTTCAACnTCACCAAGAGCTGTCTGTGGTCTTTGACCCA 



GGGACTGGAGCGTGGTCTCCTTATGATAATTGAAGGACAGGTTGAGTCTGCGC 
GCGACTGGAGCGTGGTCTTCTTATGATAATTGAAGGACAGGTTGAGTCTGCGC 

GAGCCGGGGGTGCCCAGTGTCACAGTCTGCCTCTACCTCCCCTGTGGCAGCCGC 
GAGCCGGGGGTGCCCAGTGGCACAGTCTGCCTCTACCTCCCCTGTGGCAGCCGC 

CGGCCAGGGCAAAAAAGCTCACGAAGCCGATGCTG1TGCCGCTGAGGTCCAG 
CGGCCAGGGCAAAAAAGCCCACGAAGCCGATGCTGTTGCCGCTGAGGTCCAG 

AGCCGCACATAGCGGGAGCCTGTFGGGCCAGCAGGAAGCTGGCACGCAGGAG 
AGCCGCACATAGCGGGAGTCTGTTGGGCCAGCAGGAAGCTGGCACGCAGGAG 



The polymorphic base is in bold, the additional inserted 3' mismatch is in lowercase, the amplicon size in base pairs (bp) and annealing temperature [°C] used for 
the qPCT are shown below the SNP name. 
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differentiate the products according to their size (Figure 2) 
and to genotype each sample. Samples were analysed using 
the QIAxcel ScreenGel Software (ver. 1.2.0; QIAGEN) and 
the standard DM500 analysis protocol. 

For one TLR4 SNP (T23C) we could not successfully 
optimise a set of allele-specific primers which could suc- 
cessfully genotype this base position. Both primers with 
and without the additional artificial base pair mismatch 
were tried, but an assay which would amplify DNA from 
dogs with the C allele could not optimised. Therefore, 
a Taq-man (hydrolysis) probe-based qPCR assay was 
designed using Primer3 and optimised against the syn- 
thetic positive control and residual DNA from the se- 
quencing experiments. The assay comprised a forward 
primer (GCACAGAAAATGCCAGAATG), reverse primer 
(AGGTGAACCGACATGAAACC) and two allele-specific 
probes (FAM-CTACCCGCCTGGCTGGGATAT-BHQ-1) 
and (HEX-CTACCCGCCTGGTTGGGATAT-BHQ-1). 
The genotyping was performed using GoTaq Colourless 
Master Mix (Promega) with 200nM primers, lOOnM of 
each probe and 4.5 mM Mg in a final volume of 25 ul. 
The same thermocycling protocol was used as for the 
allele-specific genotyping with a 60°C annealing tem- 
perature, but 40 cycles were performed with genotyp- 
ing based on end-point fluorescence using the MxPro 
software. Reactions were run in duplicate for each sample 
with homozygous, heterozygous and no-template controls 
were run with each reaction. 

Association study 

Statistical analysis was performed with the whole gen- 
ome association analysis toolset PLINK (http://pngu. 
mgh.harvard.edu/~purcell/plink/). Minor allele frequen- 
cies were calculated for each SNP and case-control com- 
parisons were analysed for all breeds together (using 
contingency tables and Chi-squared tests) and separately 
(using two-tailed Fisher exact tests). Statistical significance 
was set at p <0.05. 
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